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Abstract. An empirical relation between the broad line 
region (BLR) size and optical continuum luminosity is of- 
ten adopted to estimate the BLR size and then the black 
hole mass of AGNs. However, optical luminosity may not 
be a good indicator of photoionizing luminosity for ex- 
tremely radio-loud AGNs because the jets usually con- 
tribute significantly to the optical continuum. Therefore, 
the black hole masses derived for blazar-type AGNs with 
this method are probably overestimated. Here we first 
derived a tight empirical relation between the BLR size 
and the H^ emission line luminosity, i?(light — days) = 
24.05(LH3/1042ergs s-i)"-^^, from a sample of 34 AGNs 
with the BLR size estimated with the reverberation map- 
ping technique. Then we applied this relation to estimate 
the black hole masses of some AGNs and found that for 
many extremely radio-loud AGNs the black hole masses 
obtained with the R — L^f, relation are systematically 
lower than those derived previously with the R—L •> re- 

^ 5100A 

lation, while for radio-quiet and slightly radio-loud AGNs 
the results obtained with these two methods are almost 
the same. The difference of black hole masses estimated 
with these two relations increases with the radio-loudness 
for extremely radio-loud AGNS, which is consistent with 
the fact that their equivalent widths of H^ emission line 
become smaller at higher radio-loudncss. If the small H/j 
equivalent widths of extremely radio-loud AGNs are in- 
deed caused by the beaming effect, we argue that the op- 
tical emission line luminosity may be a better tracer of 
ionizing luminosity for blazar-type AGNs and the black 
hole mass derived with the R — Lh^ relation are probably 
more accurate. 

Key words: black hole physics - galaxies: active ~ galax- 
ies: nuclei - quasars: general ~ quasars: emission lines 



1. Introduction 

Supermassive black hole is essential for AGNs activities 
(Lynden-Bell 1969; Rees 1984). The black hole masses 
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of some nearby AGNs have been recently estimated by 
the reverberation mapping technique (Wandel, Peterson 
& Malkan 1999; Ho 1999; Kaspi et al. 2000), with which 
the size of the broad line region (BLR) can be measured 
from the time delay between the fiux variations of the 
continuum and the emission lines of AGNs. The black 
hole mass is then estimated using the Virial theorem from 
the BLR size and the characteristic velocity (determined 
by the full width at half-maximum (FWHM) of emission 
line). So far, the reverberation studies have yielded the 
black hole masses of about 20 Seyfert 1 galaxies and 17 
nearby bright quasars. 

An empirical relation between the BLR size (R) and 
the optical continuum luminosity at 5100A has 
been derived by Kaspi et al. (2000) using the observed 
data of 34 nearby AGNs. Because the measurement of 
the BLR size with the reverberation mapping technique 
needs long-term monitoring of continuum and emission 
line fiuxes, it is impractical for most AGNs. Therefore, 
the empirical relation has been frequently adopted to es- 
timate the BLR size and then derive the black hole masses 
for AGNs in some samples of mostly radio-quiet objects 
(Laor 2000; McLure & Dunlop 2001; Wandel 2002), and 
of purely radio-loud objects (Lacy et al. 2001; Gu, Cao & 
Jiang 2001; Oshlack, Webster & Whiting 2002). However, 
the optical luminosity of some radio-loud AGNs (espe- 
cially blazars), may not be a good indicator of ionizing 
luminosity, which is usually related to the UV/optical ra- 
diation from the accretion disk around the central black 
hole. The relativistic jets of blazar-type AGNs not only 
dominate the radio and high energy X-ray and 7 ray radi- 
ation, but also significantly contribute to the optical lumi- 
nosity in some cases (Scarpa & Urry 2002). For example, 
many optical jets have been discovered recently in AGNs 
by the HST (Scarpa et al. 1999; Jester 2003; Parma et 
al. 2003), which clearly suggests that the jets contribute 
significantly in the optical band. Furthermore, optical syn- 
chrotron radiation has been detected for many other radio- 
loud AGNs (Whiting, Webster & Francis 2001; Chiaberge, 
Capetti, & Celloti 2002; Cheung et al. 2003). Therefore, 
the measured optical continuum luminosity of some ex- 
tremely radio-loud AGNs is significantly contributed by 
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the optical radiation from the jets and may be much larger 
than the ionizing luminosity required to produce broad 
emission lines. Using the empirical relation between the 
BLR size and optical luminosity at 5100A, which was ob- 
tained based on the sample of mostly radio-quiet AGNs 
(Kaspi et al. 2000), one would significantly overestimate 
the actual BLR size and hence the black hole mass of these 
radio-loud AGNs. Oshlack et al. (2002) have shown that 
their estimated black hole masses would be lower if the 
synchrotron contribution to the optical flux is subtracted. 
However, it is not easy to make such a correction for a 
large sample of radio- loud AGNs. In addition, the contri- 
bution of the host galaxy to the optical continuum should 
also be taken into account especially when the host galaxy 
of AGNs can be resolved optically. Therefore, optical lu- 
minosity may not be a good indicator of photoionization 
luminosity of AGNs in some cases. 

In this paper we will first derive an empirical relation 
between the BLR size and the emission line luminosity 
for 34 AGNs in the sample of reverberation mapping stud- 
ies. We then argue that the BLR size obtained from the 
luminosity is more reasonable at least for some extremely 
radio-loud AGNs. Finally we apply this new empirical re- 
lation to estimate the black hole masses of some quasars 
and compare them with previous results. 

2. The relation between the BLR size and Hp 
luminosity 

Kaspi et al. (2000) have compiled the observational data 
of 17 Seyfert galaxies (Wandel, Peterson & Malkan 1999) 
and 17 nearby quasars with black hole masses estimated 
with the reverberation mapping technique. Using a linear 
fit to the available data with errors, they got an empirical 
relation between the BLR size and the optical continuum 
luminosity at 5100A as: 



i?Bifl(liglit - days) = (32.9t2;0) 



[L « /lO^ergs s'^] 

^ 5100^' ^ ' 



-ll0.700±0.033 



(1) 



With an ordinary least square (OLS) bisector method 
(Isobe et al. 1990), we can obtain such a relation 
as i?BL«(light - days) = Sl.liL^^^^^/lO^^crgs s-i]0-™i, 
which is almost the same as that shown above. We should 
keep in mind that this R — L ° relation was obtained 

. . . 5100/1 

with mostly radio-quiet AGNs. The optical continuum lu- 
minosity may not be a good indicator for photoionizing 
luminosity for some extremely radio-loud AGNs. Instead, 
the emission line luminosity may be a better indicator be- 
cause it is free from the beaming effects of the jet. 

Using the available data of BLR sizes and fluxes 
for 34 AGNs in the reverberation mapping studies, we can 
investigate the relation between the BLR size and the 
emission line luminosity (including both broad and nar- 
row components). In Table 1 we listed the BLR size and 
the luminosity data of these 34 AGNs. The luminosity 
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Fig. 1. The relations of BLR size and luminosity for 34 

AGNs in the reverberation mapping studies. The open 
and filled symbols denote Seyfert galaxies and quasars re- 
spectively. The upper panel shows the R — Lhs relation. 



The line shows the OLS bisector fit to the data. The lower 
panel shows the R — L ° relation, which is identical to 

^ 5100A ' 

the Figure 6 in Kaspi et al. (2000) except that we include 
Mrk 279 (shown as open triangle) . The line represents the 
linear fit given in Kaspi et al. (2000). 



is calculated from the flux which is available for 16 PG 
quasars, 8 AGNs in Ohio sample, and other 9 Seyfert 1 
galaxies (see references listed in Table 1). Because there 
is no available data of Hjs flux for PG 1351-1-640 (Kaspi 
et al. 2000), we exclude this object from our investiga- 
tion. In addition, we add another Seyfert galaxy Mrk 
279 in our sample because both the BLR size and the 
Hp flux have been measured recently (Santos-Lleo et al. 
2001). All luminosity data have been corrected for 
Galactic extinction using the values from NED^ (sec also 
Burstein & Heiles 1982). The cosmology with Hubble con- 
stant Hq = 75km s~^Mpc~^ and deceleration parameter 
qo = 0.5 were adopted through out the paper. 
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With these data, we derive an empirical relation be- 
tween the BLR size and luminosity. With the OLS 
bisector method we obtained: 

Log R (light - days) = (1.381 ± 0.080) + 

(0.684 ± 0.106)Lo.g (i^^/lO^'^ ergs s^^). (2) 

The slope of this relation is slightly flatter than that of 
R—L^_^^^^ relation given in Kaspi et al. (2000), consistent 

with Lh^ oc obtained by Ho & Peng (2001) for PG 

quasars. The Spearman's rank correlation coefficient of 
our R — relation is 0.91, slightly higher than 0.83 
for the R — ^^^^^^^ relation (Kaspi et al. 2000), which 
implies that the R — Lh^ relation is slightly tighter than 
the R — L ' relation. We also used the linear fit to the 

5100/1 

data with errors that was used by Kaspi ct al. (2000) to 



derive the R — L ° relation. Applying this method to 

5100A " " ■ ° 

the R — LHf, relation yields the same result as the OLS 
bisector method. 

In Fig. 1 we show the dependence of the BLR size 
on L Ha and L - . The two relations are similar and 
thus the R — Lhi, relation can be an alternative of the 
R — L ° relation in estimating the BLR size for radio- 

5100A ° 

quiet AGNs. 



3. Comparison of black hole mass estimation of 
AGNs from two relations 

Since the R — L ° relation obtained by Kaspi et al. 

5100A 

(2000) has been frequently used to estimate the BLR size 
and the black hole mass of both radio-quiet and radio-loud 
AGNs, it is important to investigate the applicability of 
such an approach for radio-loud objects. 

Brotherton (1996) studied the emission line properties 
of 59 radio- loud quasars. We adopted his published val- 
ues of absolute V-band magnitude, equivalent width and 
FWHM of Hjs emission line. The continuum luminosity 
at 5100A and the Hp luminosity (scaled to our cosmology 
parameters) were calculated after considering the Galac- 
tic extinction and K-correction (optical spectral index was 
assumed to be 0.3). We then estimated the BLR size using 
both R — Lh„ and R — L ° relations and derived the 

'3 5W0A 

black hole mass with the formula Mbh = 3V^^^^i?/4G' 
(here we assumed the BLR velocity V ^ {V5/2)Vfwhm 
as in Kaspi et al. 2000). With these two relations, we also 
estimated the black hole masses of another 27 radio-loud 
quasars with available data of both the equivalent width 
and FWHM of Hf3 emission line in the Parkes Half- Jansky 
flat-spectrum Sample (PHFS) (Drinkwater et al. 1997; 
Francis, Whiting & Webster 2000; Oshlack ct al. 2001). 
We compared the black hole masses obtained with the 
R — LHfi and R — L^,^^^^ relations in Figure 2. Evidently 
the masses obtained with the R — L^f,^ relation are sys- 
tematically lower that those obtained with the R — L 
relation for some extremely radio-loud quasars. 
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Fig. 2. Upper panel: Comparison of the black hole masses 

of radio-loud quasars estimated with R — Lh^^ and R — 
"^siool relations. The squares represent 59 quasars in 
Brotherton (1996) and the triangles represent 27 PHFS 
quasars in Oshlack ct al. (2001). The diagonal line shows 
the relation where both masses are identical. Lower panel; 
The ratios of black hole masses estimated with two dif- 
ferent relations are plotted against the radio-loudncss of 
radio-loud quasars. The dashed line indicates the case 
where the two black hole mass estimates are identical. Ev- 
idently MBjy[5100A] becomes systematically larger than 
MbhIH^}] at higher radio-loudness. 



In Figure 2 we also show how the difference of black 
hole masses obtained with these two relations varies with 
the radio-loudness for these two samples of radio-loud 
quasars. The values of radio-loudness of 27 PHFS quasars 
were taken from Oshlack et al. (2001) and those of 59 
quasars in Brotherton (1996) were calculated from the Ry 
value (defined as the ratio of core radio luminosity and the 
V-band optical luminosity) and the core and extended ra- 
dio luminosity values listed in his Table 1. It is clear that 
the difference of black hole masses is small when the radio- 
loudness is small but becomes larger as the radio-loudness 
increases. The Spearman's rank correlation coefficient is - 
0.34, implying a modest correlation between the difference 
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Fig. 3. The relation between the equivalent width of 
emission line and the radio-loudness of radio-loud quasars. 
The symbols have the same meaning as in Figure 2. 



of black hole masses and the radio-loudness. For some in- 
dividual quasars with higher radio-loudness, the black hole 
mass estimated with the R— L ° relation can be 3~10 

5100A 

times larger than that estimated with the R—Lh^ relation. 
In Figure 3 we also plotted the equivalent width (EW) of 

emission line against the radio-loudness for objects in 
these two radio-loud AGN samples. We can see that the 
EW(H/3) becomes smaller at higher radio-loudness. Such 
a relation, although with only a modest Spearman's rank 
correlation coefficient of -0.33, indicates that the smaller 
EW(H/3) of some extremely radio-loud AGNs could be at 
least partly due to the beaming effects. 

For radio-quiet AGNs, however, both the optical con- 
tinuum and emission line luminosities are probably free 
from the jet contributions and therefore both can be good 
tracers of photoionization luminosity. We check this by 
using the data of 70 low-redshift radio-quiet quasars in 
the Palomar-Green survey. The emission line properties 
of these quasars have been studied by Boroson & Green 
(1992). The contimmm luminosity at 5100^1 and the Hjj 
luminosity were estimated from the absolute V-band mag- 
nitude [Mv) and the equivalent width of Hp emission line 
listed in their Table 1 and Table 2. Using the R — Lh.^ and 
R—L " relations we estimated the BLR sizes and black 

5W0A 

hole masses of these radio-quiet quasars. The results from 

the two relations are almost identical (see Figure 4). This 
is also indicated by the normalized -y^ value of the devi- 
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Fig. 4. Upper panel: Comparison of black hole masses 
of 70 radio-quiet PG quasars estimated with the R — 
and R — L ° relations. The diagonal line shows the case 

5100A ° 

where the masses are identical. Lower panel: The ratios 
of black hole masses estimated with two different rela- 
tions are plotted against the radio-loudness of radio-quiet 
quasars. The dashed line indicates that the case where the 
two black hole mass estimates are identical. 



ation of the points plotted in the upper panel of Figure 
4, which is 0.94, much smaller than the value 2.58 for the 
points plotted in the upper panel of Figure 2 for radio-loud 
AGNs. From the lower panel of Figure 4, we can also see 
that the difference of the two black hole mass estimates 
does not correlated with the radio-loudness for radio-quiet 
quasars. The Spearman's rank test gives a correlation coef- 
ficient of only 0.05, much smaller than that for radio-loud 
AGNs. 

4. Discussions 

Using the empirical relation between the BLR size and 

optical continuum luminosity possibly induces an overes- 
timation of the BLR size and hence the black hole mass of 
some extremely radio-loud AGNs because of the jet con- 
tribution to the optical luminosity. We derived another 
empirical relation between the BLR size and the emission 
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line luminosity, and demonstrated that it can be used to 
estimate the BLR size and black hole mass of both radio- 
quiet and radio-loud AGNs. If the relativistic jets and host 
galaxy have significant contributions to the optical con- 
tinuum, the emission line luminosity is probably a bet- 
ter tracer of ionizing luminosity. Comparisons of the esti- 
mated black hole masses with these two different empirical 
relations clearly indicate that the difference becomes sig- 
nificant if the radio-loudness of AGNs is larger. Using the 
R — Lh,^ relation may result in more accurate estimations 
of black hole masses of some blazar-type AGNs. 

In our study we focused on the possible effects of beam- 
ing on the optical continuum in radio-loud AGNs and ig- 
nored the difference in BLR physics between radio-loud 
and radio-quiet AGNs. The modest correlation between 
the EW(H/3) and the radio-loudness may indicate the pres- 
ence of beaming effects, though some other effects such as 
a lower covering factor of the BLR of radio-loud AGNs can 
also lead to smaller EW (H^^) values. Because currently we 
know little about the difference of BLR physics between 
radio-loud and radio-quiet AGNs, to prove the validity 
of our approach it is necessary to compare our estimated 
black hole mass with an independent estimate, for exam- 
ple, from the correlations of black hole mass with central 
velocity dispersion and host galaxy luminosity. Unfortu- 
nately, not many measured values of central velocity dis- 
persion or host galaxy luminosity for extremely radio-loud 
AGNs are available. Although there are 8 objects in the 
sample of Brotherton (1996) with measured host magni- 
tude (McLure & Dunlop 2001), the radio-loudness of these 
objects are mostly smaller than 1000 and thus the differ- 
ence estimated with the R — Lh„ and R — L ° relations 

1^ 5100A 

is rather small. The velocity dispersion measurements for 

radio-loud quasars are not available and the [OIII] pro- 
file in radio-loud AGNs may not be adopted to estimate 
the central velocity dispersion because of its complexity. 
Therefore, further imaging studies on the host galaxy and 
spectroscopic measurements of the central velocity disper- 
sions of a large sample of extremely radio-loud quasars are 
still desired to confirm our results. 

The advantage of using the R — Lh^ relation is that 
we can estimate the black hole mass of AGNs with only 
two observed parameters, namely the H^^; line luminosity 
and its FWHM, and it can be applied to a larger sample 
of AGNs with redshift smaller than 0.8. In principle, one 
can analogously investigate the relation between the BLR 
size and the luminosity of some ultraviolet emission lines 
such as Mgll and CIV, which may be used to estimate the 
black hole mass of some high redshift AGNs. Some recent 
studies have suggested to use the ultraviolet continuum 
luminosity and the FWHM of ultraviolet emission lines 
to estimate the black hole mass of high redshift AGNs 
(Vestergaard 2002; McLure & Jarvis 2002). However, the 
ultraviolet continuum luminosity can similarly suffer the 
serious contaminations from jet and Blamer continuum, 
therefore the luminosity of ultraviolet emission line again 



may be a better indicator of ionizing luminosity than the 
ultraviolet continuum luminosity. 

Finally, one should be cautious to the uncertainties in 
estimating the black hole mass of AGNs using the R — Lh^ 
relation. Firstly, the variations of H^ emission line flux and 
its FWHM are common in AGNs. Estimating the black 
hole mass with the values of these two parameters in a 
single spectrum may lead to large errors. Secondly, the 
different inclination of the BLR may also significantly af- 
fect the results (McLure & Dimlop 2001; Wu & Han 2001). 
If the BLR has a flatten geometry and the inclination of 
BLR is rather small, our derived values of black hole mass 
may significantly underestimated. However, the ratio of 
black hole masses estimated with two empirical relations 
does not depend on the inclination. Better understandings 
of BLR geometry and dynamics arc absolutely needed to 
diminish the uncertainties in deriving the black hole mass 
of AGNs (Krolik 2001). 
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Name 


z 


Ab 


Rblr 
(light-days) 


L 

.5100A 

(10*^erg/s) 


(10*^ erg/s) 


Ref 


3C 120 


0.033 


0.570 


42 0+^^-" 


73.00 


± 


13.00 


1.222 


± 


0.120 


1 


3C 390.3 


0.056 


0.170 


00 q+6-3 


64.00 


± 


11.00 


1.433 


± 


0.171 


3 


Akn 120 


0.033 


0.400 


37.4tl:l 


139.00 


± 


26.00 


2.496 


± 


0.463 


1 


F 9 


0.046 


0.000 


16.3+?-^ 


137.00 


± 


15.00 


2.056 


± 


0.174 


4 


IC 4329A 


0.016 


0.000 


lA+li 


16.40 


± 


2.10 


0.183 


± 


0.014 


5 


Mrk 110 


0.035 


0.000 


18.81^;^ 


38.00 


± 


13.00 


0.811 


± 


0.203 


1 


Mrk 279 


0.030 


0.000 


16.21^-1 


66.00 


± 


6.30 


1.023 


± 


0.071 


12 


Mrk 335 


0.026 


0.100 




62.20 


± 


5.70 


1.160 


± 


0.073 


1 


Mrk 509 


0.034 


0.180 


76.7tt:l 


147.00 


± 


15.00 


3.070 


± 


0.332 


1 


Mrk 590 


0.026 


0.050 




51.00 


± 


9.60 


0.498 


± 


0.143 


1 


Mrk 79 


0.022 


0.230 




42.30 


± 


5.60 


0.619 


± 


0.041 


1 


Mrk 817 


0.031 


0.000 




52.60 


± 


7.70 


0.740 


± 


0.123 


1 


NGC 3227 


0.004 


0.020 


10 Q+5-6 

lU.a_j^Q g 


2.02 


± 


0.11 


0.017 


± 


0.002 


6 


NGC 3783 


0.010 


0.470 


4.5ilf 


17.70 


± 


1.50 


0.292 


± 


0.021 


7 


NGC 4051 


0.002 


0.000 


6.5+ti 


0.525 


± 


0.03 


0.0048 


± 


0.0005 


8 


NGC 4151 


0.003 


0.000 


s.o+li 


7.20 


± 


0.42 


0.152 


± 


0.010 


9 


NGC 5548 


0.017 


0.000 


21.2+tv 


27.00 


± 


5.30 


0.421 


± 


0.092 


10 


NGC 7469 


0.016 


0.120 


4.9+?f 


55.30 


± 


1.60 


0.423 


± 


0.019 


11 


PG 0026 


0.142 


0.130 


113.01^?;° 


700.00 


± 


100.00 


5.693 


± 


0.493 


2 


PC 0052 


0.155 


0.120 


134.01^1;° 


650.00 


± 


110.00 


9.595 


± 


1.119 


2 


PG 0804 


0.100 


0.110 


I56.0t\l° 


660.00 


± 


120.00 


13.95 


± 


0.908 


2 


PG 0844 


0.064 


0.080 


24.2iiV 


172.00 


± 


17.00 


2.585 


± 


0.257 


2 


PG 0953 


0.239 


0.000 




1190.00 


± 


160.00 


19.39 


± 


1.129 


2 


PG 1211 


0.085 


0.130 


lOl.Olg:" 


493.00 


± 


80.00 


8.588 


± 


1.056 


2 


PG 1226 


0.158 


0.000 


387.o1^«:E; 


6440.00 


± 


770.00 


88.37 


± 


7.192 


2 


PG 1229 


0.064 


0.000 




94.00 


± 


10.0 


1.601 


± 


0.202 


2 


PG 1307 


0.155 


0.020 


124.01^^0 


527.00 


± 


52.00 


9.603 


± 


1.301 


2 


PG 1411 


0.089 


0.000 


102.0l^«;0 


325.00 


± 


28.00 


5.268 


± 


0.285 


2 


PG 1426 


0.086 


0.120 




409.00 


± 


63.00 


4.952 


± 


0.465 


2 


PG 1613 


0.129 


0.040 


39.0l?«o 


696.00 


± 


87.00 


7.014 


± 


0.451 


2 


PG 1617 


0.114 


0.150 




237.00 


± 


41.00 


4.060 


± 


0.493 


2 


PG 1700 


0.292 


0.020 


oo r|+190.0 


2710.00 


± 


190.00 


35.29 


± 


1.851 


2 


PG 1704 


0.371 


0.000 




3560.00 


± 


520.00 


9.752 


± 


1.257 


2 


PG 2130 


0.061 


0.170 


200.0+%° 


216.00 


± 


20.00 


4.241 


± 


0.381 


2 



Notes: The Galactic extinction values are adopted from NED. Data of Rblr and L ° are taken from Kaspi et al. (2000). 
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The luminosity is calculated from the flux given in the literature (sec the column "ref"). References: (1) Peterson et al. 
1998; (2) Kaspi et al. 2000; (3) Dietrich et al. (1998); (4) Santos-Lleo et al. (1997); (5) Winge et al. (1996); (6) Winge et al. 
(1995); (7) Strirpe et al. (1994); (8) Peterson et al. (2000); (9) Kaspi et al. (1996); (10) Peterson et al. (2002); (11) Collier 
et al. (1996); (12) Santos-Lleo et al. (2001). 



